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Abstract. We describe the results of recent studies of a chromomagnetic instability and a gluonic 
phase in neutral two-flavor quark matter at nonzero temperature. 
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INTRODUCTION 

Sufficiently cold and dense quark matter is a color superconductor. Color superconduct- 
ing phases at "moderate" density, which is relevant for quark matter in the interior of 
compact stars, have a rich phase structure and have recently attracted a great deal of 
interest yj]. 

Bulk matter in the interior of compact stars should be in /3 -equilibrium and be neutral 
with respect to electric and color charges. In the two-flavor case, these conditions sepa- 
rate the Fermi momenta of up and down quarks and, as a consequence, the ordinary BCS 
state (2SC) is not always energetically favored over other unconventional states, e.g., the 
gapless two-flavor color superconductivity (g2SC) Q2Q. However, the 2SC/g2SC phases 
suffer from a chromomagnetic instability, indicated by imaginary Meissner masses of 
some gluons yfl. The instability related to gluons of color 4-7 occurs when the ratio 
of the gap over the chemical potential mismatch, A/ 5/1, decreases below a value y2. 
Resolving the chromomagnetic instability and clarifying the nature of the true ground 
state of dense quark matter are central issues in the study of color superconductivity 

i,a a a s, §, y, o, gum h h em qi gi m . 



So far, several candidates for the ground state have been proposed JUS 6 ,0,0, Gil]- 
In this work, we study the gluonic phase (gluonic cylindrical phase II) [9] at nonzero 
temperature. 

FREE ENERGY OF GLUONIC PHASE 
Model 

In order to study the gluonic phase, we use a gauged Nambu-Jona-Lasinio model with 
massless up and down quarks: 

££ = yW+fLfty+GD ^ir 5 £e b Cxfr T ^j (yCiy 5 ££ b y) - l -F^ v F a ^\ (1) 



where the quark field \j/ carries flavor (i,j= l,---,Nf with Nf = 2) and color (oc,/3 = 
l,---,N c with N c = 3) indices, C is the charge conjugation matrix; (e)' k = £ lk and 
_ £ bafi ^ ^ antisymmetric tensors in flavor and color spaces, respectively. 
The diquark coupling strength in the scalar (J p = + ) color- antitriplet channel is denoted 
by Go- The co variant derivative and the field strength tensor are defined as 

Dp = dfi — igA^T", F« v = d^A a v - d v Al +g f^A b fl A c v . (2) 

The elements of the diagonal matrix of quark chemical potentials (1 in /3 -equilibrated 
neutral two-flavor quark matter are given by 

Hur = Vug = /i-5jU, \L dr = ji dg = A + 5/1, 

^ = £-5/1-^8, ji db = ^ + 8ji-ji % , (3) 



with 



p = M -5t + »« M = * (4) 



In Nambu-Gor'kov space, the inverse full quark propagator S 1 (p) is written as 

(So 



s-\p)=[ <*r *w ), (5) 



with 



= f Pll + (a - 5MT 3 )y° +g^A« r , (6) 
(SoT 1 = tvia - (A - s^V -gy^r r , (7) 



and 



4> = -z'££°y 5 A, 4> + = -iee b y 5 A.. (8) 



Here T 3 = diag(l, —1) is a matrix in flavor space. 

For the gluonic phase, B = (gA^) is the most relevant vector condensate, because the 
tachyonic mode related to gluons 4-7 exists in the direction of B MAIM- Besides B, to 
ensure color neutrality at B ^ we have to introduce /I3 = (gA^) [9]. In what follows, 
however, we neglect color chemical potentials /ig and /I3, since their effect on the free 
energy is negligibly small for a s ~ 1 . Taking account of B, the free energy of the gluonic 
phase in the one-loop approximation is given by 



_^__ly f d 3 p 
+ 4G D 2^.1 {2k) 



(9) 



where /3 = l/T , the £ a 's are quasi-quark energies and the sum runs over all particle and 
anti-particle £ a 's. Note that the e a 's depend on B through the co variant derivatives in the 
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FIGURE 1. Left: The Ao dependence of the gap A, the chemical potential mismatch 5/x, and the vector 
condensate B at T = 0. Right: The free energy (measured with respect to the 2SC/g2SC/NQ phases) as a 
function of B at T = 0. The results are plotted for = 400 MeV. 



quark propagator. Here, we added a contribution from electrons (first line on the r.h.s.). 
In order to evaluate loop diagrams we use a three-momentum cutoff A = 653.3 MeV 
throughout this work. In addition, in order to remove a cutoff dependence of the free 
energy we introduce the following subtraction: 

Q R = Q(A, 8n,B;[i,T) -£2(0,0,5; 0,0). (10) 

Note that this free energy subtraction is not adequate to remove the cutoff dependence of 
the free energy at T > [fl5ll . However, we have carefully checked that it is nothing but 
a cutoff artifact and moreover is negligibly small at /i ~ 400 MeV and at temperatures 
of interest (20 MeV at most). In order to find the gluonic phase, we first solve a set of 
coupled equations, 

as a function of B and, then, compute Q,r(B). Finally, the minimum of Q.r(B) determines 
the gluonic phase. 



Numerical results 

In the left panel of Fig. \T\ we plot A, 5/i and B as a function of Ao, the 2SC gap 
at 5/1 = 0. (Essentially, Ao corresponds to the diquark coupling strength.) We find that 
the gluonic phase exists in the window 66 MeV < Ao < 162 MeV and that a strongly 
(weakly) first-order transition between the gluonic phase and the NQ (2SC) phase takes 
place at Aq = 66 (162) MeV. Note that the results are consistent with those in Ref. 
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FIGURE 2. The temperature dependence of the free energy as a function of B for Ao = 75 MeV (left) 
and Aq = 85 MeV (right). The quark chemical potential is taken to be jA — 400 MeV. 



UM, where the color chemical potentials /I3 and were treated self-consistently. 
The right panel of Fig. \T\ shows the free energy Q.r(B) (measured with respect to the 
2SC/g2SC/NQ phases at B = 0) as a function of B at T = 0. One can see that, in the 
cases of Ao = 75, 85 MeV, the gluonic phase is energetically more favored than the 
chromomagnetically stable NQ phase at B = 0. [For small values of B, we found that the 
system is in the ungapped phase and the free energy behaves like £Ir(B) ~ 0(B 4 ).] On 
the other hand, for A = 100 MeV and A = 140 MeV, the 2SC/g2SC phases at B = 
are unstable and, then, the gluonic phase is energetically favored, as expected. 

In Fig. [2l we plot the temperature dependence of the free energy for Ao = 75 MeV 
(left) and Ao = 85 MeV (right). At Ao = 75 MeV, the free energy gain gets rapidly 
reduced when T is increased, but the vacuum expectation value of B is only slightly 
reduced. Then we find a strongly first-order transition at T ~ 14 MeV. In the case of Ao = 
85 MeV, the gluonic phase is energetically more favored than the chromomagnetically 
stable NQ phase at low temperature. At T ~ 9 MeV, the stable NQ phase undergoes a 
phase transition into the unstable g2SC phase. The gluonic phase is energetically favored 
over the unstable g2SC phases until the temperature reaches T ~ 20 MeV. Above this 
temperature, the gapless dispersion relation is smoothed out and the stable g2SC phase 
is favored. 

Here, we would like to make a comment regarding the order of the phase transitions. 
In the weak coupling regime, the transition from the gluonic phase to the NQ phase 
is strongly of first order. On the other hand, the phase transitions from the gluonic 
phase to the 2SC/g2SC phases in the intermediate coupling regime are likely to be of 
second order. However, evaluating the free energy near the critical temperatures self- 
consistently is not easy. Thus we cannot exclude the possibility of weakly first-order 
transitions. 
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FIGURE 3. Left: The Ao-T phase diagram of electrically neutral two-flavor quark matter. At T = 0, 
the g2SC phase phase exists in the window 92 MeV < Ao < 134 MeV and the 2SC window is given by 
Ao > 134 MeV. The unstable region, where gluons 4-7 are tachyonic, is depicted by the region enclosed 
by the thick solid line. The g2SC phase and a part of the 2SC phase (92 MeV < A < 162 MeV) suffer 
from the chromomagnetic instability at T = 0. Right: Schematic phase diagram of the gluonic phase ll20ll . 
In the region enclosed by the thick solid and dotted lines, the gluonic phase is more favored than the 
2SC/g2SC/NQ phases. First-order phase boundary is indicated by the thick dotted line. The thick solid 
line denotes the line of second-order or weakly first-order transitions. 



CONCLUSION 

In the left panel of Fig.[3l we plot the phase diagram of the electrically neutral 2SC/g2SC 
phases in the plane of Ao versus T. The region enclosed by the thick solid line is unstable 
because the Meissner mass of gluons 4-7 is imaginary there. Note, however, that we 
did not examine the global structure of the free energy self-consistently. Therefore we 
should regard this result as the tendency toward the gluonic phase in the 2SC/g2SC 
phases. In fact, at T = 0, the gluonic cylindrical phase II exists in the wider window 
66 MeV < A < 162 MeV (see also Ref. (Hi). 

In this work, we studied the phase structure of the gluonic cylindrical phase II at 
nonzero temperature. In the weak coupling regime, we found that the phase transition 
between the gluonic phase and the NQ phase is of strongly first order. In contrast, in the 
intermediate and strong coupling regimes, the transition is of second order or weakly 
first order. Therefore, in these regimes, the critical line shown in the left panel of Fig. [3] 
is not altered much by the self-consistent analysis. In the right panel of Fig. [31 we show 
the schematic phase diagram of the gluonic phase. One can see that the gluonic phase 
wins against a part of the NQ phase and enlarges its region. The phase transition to the 
NQ phase is of strongly first order (dotted line) while the transition to the 2SC/g2SC 
phases is of second order or weakly first order (solid line). The critical temperature 
for the gluonic phase is roughly given by T c ~ 20 MeV. (It is known that the gluonic 
phase survives until the temperature reaches T c ~ where 5jU^ N ^ denotes 

the chemical potential mismatch in the NQ phase at T = lfl5|l . The present model 



parameters yield 5jU (NQ) ~ 40 MeV and, hence, T c ~ 20 MeV.) 

Finally, we comment on the outlook for future studies. Several candidates for the true 
ground state (for instance, other types of the gluonic phase [9|] and a LOFF state with 
multiple plane waves fl]) have been proposed. It is then quite interesting to study the 
free energy of these phases and their stability against the chromomagnetic instability. In 
addition, since the gluonic phase dominates the low-temperature region of the phase di- 
agram for the wide range of the coupling strength, it would have important implications 
for the properties of compact stars. 
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